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Abstract. 

We investigate the roughening of shear cracks running along the interface between 
a thin film and a rigid substrate. We demonstrate that short-range correlated 
fluctuations of the interface strength lead to self-affine roughening of the crack front as 
the driving force (the applied shear stress/stress intensity factor) increases towards 
a critical value. We investigate the disorder-induced perturbations of the crack 
displacement field and crack energy, and use the results to determine the crack pinning 
force and to assess the shape of the critical crack. The analytical arguments are 
validated by comparison with simulations of interface cracking. 



1. Introduction 

Strength of materials is governed by the properties of defects. In particular, in order 
to understand materials failure, one needs to study the nucleation of new and the 
propagation of existing cracks. While many investigations (for an overview, see PQ) 
have focused on cracks in bulk materials, in the present work we consider near-surface 
cracks running along a weak interface which separates a thin elastic film from a rigid 
substrate and acts as a prescribed fracture surface. Failure is induced by shear tractions 
that are applied to the free surface of the film and induce shear stresses in the film and 
on the interface. The considered geometry applies to a variety of materials problems, 
such as shear failure and/or abrasive wear of coatings and shear- induced delamination 
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of thin films. One may also consider geophysical applications where the load acting on 
the weak interface originates from gravitational body forces, such as in the initiation of 
snow slab avalanches by failure of weak interfaces in snow stratifications j2j [3l HI [5] . 

While in a companion paper [6] we have addressed the problem of crack nucleation, 
in the present study we focus on the behavior of pre-existing interface cracks. We 
assume a disordered interface which exhibits a random pattern of weak and strong 
regions. Hence, the shape of an advancing crack becomes distorted as the crack advances 
across weak regions and gets arrested at regions of elevated strength. At the same time, 
elastic forces resist the distortion and try to keep the crack front straight. The interplay 
between crack front elasticity and disorder gives rise to the mutually related phenomena 
of crack front roughening and crack pinning. 

The roughening of bulk cracks has been investigated extensively in view of 
explaining the emergence of self-affine fracture surfaces (see e.g. [7j, [H]). The 
long-range elastic self-interaction of the two-dimensionally perturbed front of a bulk 
crack was analyzed by Gao and Rice [9] and the implications on quasi-static crack 
propagation in disordered materials were discussed by Ramananthan et. al. [10J. In 
this work, cracks were envisaged as one- dimensional elastic manifolds with long-range 
self-interactions, which are driven by an external force (the applied stress intensity 
factor) through a 'landscape' of short-range correlated disorder. From this point of 
view, the in-plane roughening of the crack front is expected to be governed by the same 
roughness, correlation length and dynamic exponents which characterize the depinning 
of a contact line. These exponents were evaluated by Ertas and Kardar using a one- 
loop renormalization scheme [11] . A two- loop renormalization calculation of the same 
exponents was done by Chauve et. al. [321 [13]. While results of the one-loop calculation 
yield roughness exponents that are generally too low when compared with experimental 
data, the two-loop correction leads to an in-plane roughness exponent of about ( = 0.5 
which is in line with available experimental observations. 

From an experimental point of view, the out-of-plane roughness of bulk cracks 
(i.e., the roughness of the fracture surface) has been studied in many papers, whereas 
comparatively few studies have been devoted to investigating in-plane roughness, i.e. 
the shape of the crack front in the direction of crack propagation. The obvious reason is 
that observing the crack front requires tricky in situ methods, whereas the out-of-plane 
roughness of the crack transfers to the fracture surface which can be easily accessed after 
the specimen has failed. In-plane roughening of crack fronts was studied by Delaplace et. 
al. [Hj for cracks running along the interface between two PMMA plates. These authors 
report typical roughness exponents of ( ~ 0.5 — 0.6. The burst-like dynamics of crack 
propagation was studied by Vellinga et. al. [15] who used image correlation methods to 
monitor the propagation of a crack along the interface between a glass subtrate and a 
polymer film. 

The above mentioned theoretical investigations refer to bulk cracks. In thin film 
geometries, on the other hand, the range of elastic interactions is limited by the film 
thickness and this modifies the energy functional for shear cracks [16] as well as for 
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tensile cracks [T7] and anticracks [U E] • Crack front roughening in thin film geometries 
was investigated by numerical simulation of a tensile crack running between two elastic 
sheets [32], and by analytical arguments and numerical simulation of a shear crack 
underneath a thin elastic substrate [18] . In the present paper we further develop this 
investigation and show that the roughening of crack fronts in thin film geometries is 
governed by a complicated interplay of length scales. In Section 2 we determine the 
elastic energy functional of a thin film with a general interfacial displacement field and 
use this to evaluate displacement fields and elastic energies of straight and perturbed 
mode-II cracks. We then use in Section 3 standard scaling arguments from pinning 
theory to determine the pinning length above which the crack roughens, and evaluate 
the associated increase in the critical crack driving force. Furthermore we show that 
the presence of the opposite crack front leads to emergence of a second critical length 
above which the crack is flat again. The scaling arguments are validated by comparison 
with simulation results (Section 4). 

Most investigations of crack front roughening envisage the dynamics of the crack 
front as a function of the driving force (the stress intensity factor). However, the stress 
intensity factor is itself an increasing function of the crack length and the acting stress. 
In experimental investigations the driving force can be controlled in such a manner that 
stable crack propagation takes place [HI [15]. In practice, however, systems are often 
subject to fixed loads. In this case, the increase of driving force with increasing crack 
length leads to catastrophic failure by unstable crack propagation once a critical stress 
threshold is reached. In our discussion section we assess the influence of disorder on the 
failure stress of a cracked interface under conditions of load control. We demonstrate 
that the advance of the crack front during roughening leads to complex behavior where 
increasing the disorder may either result in weakening by disorder-induced crack advance 
or in strengthening by crack pinning, while at very high disorder a change in failure mode 
is observed. 

2. Formulation of the elastic problem 

2.1. General considerations 

We consider an elastic film of thickness D tethered to a rigid substrate. The interface 
with the substrate is modeled as a cohesive layer in the plane z = 0. The response of the 
interface to shear loads is characterized by the scalar stress-displacement relationship 
o~xz(x,y,z = 0) = r(u(x,y)) where a xz is the shear stress at the interface, r is the 
interfacial shear strength, and u(x, y) = w x (x, y, z — 0) is the shear displacement across 
the interface. The maximum stress that can be supported by the interface is denoted 
as tm, and the specific fracture energy per unit interface area is given by the integral 



where u denotes the characteristic displacement-to-failure. Structural disorder of the 
interface is modelled in terms of short-range correlated fluctuations of the fracture energy 




(1) 
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Wf = Wf(x, y). The stochastic process Wf(x, y) is assumed to have finite first and second 
moments (Wf) and (W 2 ), and finite correlation length £. 

The film is loaded by spatially homogeneous tractions applied to its free surface at 
z = D, giving rise to a space-independent 'external' shear stresss cr^ T . This external 
stress superimposes on the internal stresses that are associated with the interfacial 
displacement field u(x,y). To evaluate these stresses we start from the elastic energy 
functional associated with a generic displacement vector field w(r). This functional is 
given by 

/ \ 2 
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where the Wij are the components of the strain tensor, W{j = (diWj + djWi)/2, and indices 
assume the values {i, j} G {x, y, z}. The parameter A can be expressed in terms of the 
shear modulus \x and Poisson's ratio v as A = 2[ivj (1 — 2v). Energy minimization gives 
the equilibrium equation 

V 2 w + - — ^grad (div w) = . (3) 

We have to solve this equation for the boundary conditions w x (x,y,z = 0) = 
u(x, y), w y = w z = 0. In Fourier space the solution is given by 
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where u(k x , k y ) is the Fourier transform of u(x,y) and rj — 1/(3 — Au). In the present 
work we are interested in the behavior of interface cracks underneath thin films. We 
therefore assume that the film thickness D is much smaller than the characteristic 
length of variations in the displacement field, such that \kz\ « 1 for all z < D. The 
displacement field is then approximately given by 
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Inserting this into Eqn. ([2]) and using again that kD 1 gives 
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(2 — 2u)/{l — 2u). Reverting to spatial coordinates, we find the energy 
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The total energy of the system is obtained by adding to the elastic energy H the work 
done by the shear stresses at the interface: 

ru(x) 



G(u) = ~Jdx I " (o-*r - r(u)) du 



(8) 



^ttt + IviTo + <*ZT ~ <u) = , (9) 



Minimizing the total energy functional E(u) = G(u) + H(u) leads to the equilibrium 
condition 

d 2 u j d 2 u 

' dx 2 v dy 2 

where the mode-II and mode-Ill interaction constants I x and I y are given by I x = 
Dfi(l + a) and I y = 2D\i. 

2.2. Displacement field of a mode-II shear crack and crack energy functional for 
perturbations of the crack front 

To analyse the behavior of a crack in a disordered medium, we first evaluate the 
properties of a straight mode-II crack and then study crack front roughening in terms 
of perturbations of this reference crack. 

We consider a crack of width 21 with crack fronts located at x = ±Z. In the brittle 
limit t{u) = W{S(u) (S(u) is Dirac's delta function) the corresponding displacement field 
uo is given by 

( (72 _ 2\ EXT 

Si X II** x (z [_/ 71 

u (x) = \ 2I X ' E[ ' tJ ' (10) 

[ else . 

Inserting the displacement field, Eqn. (flQl) . into the energy functional, Eq. <£7§, gives 
the elastic energy of the reference crack: 

The derivative of the elastic energy with respect to I is the crack driving force per unit 
length (the stress intensity factor). Since the crack has two fronts, the driving force for 
each front is 

(la EXT ) 2 

f(l) = [ ~^ L ■ (12) 

This driving force is resisted by the fracture energy of the interface. For the reference 
crack with fracture energy (Wf), equating / and (Wf) defines a critical crack length 

l c = sj2(Wi)I x /o^ . (13) 

We now consider situations where the fracture energy exhibits random spatial 
fluctuations. In this case, parts of the crack front may advance across regions of reduced 
fracture energy and then get trapped at locations of elevated fracture energy. Thus, the 
crack may advance and roughen before it gets critical. To study the energetics of crack 
front roughening, we envisage a generic perturbation Sl(y) of the straight crack front at 
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x = I. This corresponds to a perturbation Su(x, y) of the displacement field which along 
the front of the unperturbed crack must fulfil the boundary condition 
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where Uq is given by Eqn. (TTUj) . Over the width of the crack, the displacement field 
fulfils the partial differential equation 
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By inserting this solution into the elastic energy functional, Eqn. (171), we can evaluate 
the energy cost associated with the perturbation. In Fourier space the elastic energy 
variation associated with the crack front perturbation 5l(k y ) is given by 

2 



Ixly | k y 
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(17) 



that is, the front of a perturbed interface crack behaves like an elastic manifold with 
long-range elasticity. Similar behavior can be found for bulk cracks or contact lines 
IP- 

In deriving Eqn. ([17]) we have assumed that k y l ^> 1 such that terms of order 
can be neglected. This implies that our treatment applies only to fluctuations 
with wavelengths that are short in comparison with the crack width. In the opposite 
case y < 1, the boundary conditions at the opposite crack front (x = —I) need to be 
considered as well. Assuming that front unperturbed, the crack energy change can be 
obtained by simply inserting the modified crack length, I — > I + 5l(y), into Eqn. fill I) 
and integrating over y. For perturbations of zero average this gives in Fourier space 

m * wA^y i w^'ito • (i8) 

This implies that crack elasticity on scales above the crack width becomes of infinite 
range: the interaction kernel is independent of the wavevector. 



3. Crack pinning and crack roughening 

We now study the changes in crack geometry caused by the disorder. The characteristic 
fracture energy fluctuation experienced by a smooth crack front segment of length A is 
of the order of 

1/2 

<Wf = ( t I (Wf)C v , (19) 
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where Cy denotes the coefficient of variation of the fracture energy distribution, C v = 
(Wf)/(Wf) 2 . The correlation length of the fracture energy fluctuations perpendicular 
to the crack is given by £. Crack front roughening on scale A occurs if the energy gain 
by adjusting to the fluctuating energy profile is larger than the energy that has to be 
expended to change the crack shape. Combining Eqs. (|T9l) and (fTTl) . and considering 51 
in Eq. ffTTl) as a perturbation of amplitude £ and characteristic wavelength A, we find 
the condition 



5H{i, A) » (l^D 2 ]]^/ < £A(W f (A) = (£ 3 A) 1/2 (m>CV . (20) 
Hence, roughening occurs above a critical scale 

, ,. ('^ T )V 'y(')U (21) 



U*(W t )*<% I x 

Above this scale, the crack adjusts to the disordered pinning 'landscape'. This leads 
to an incresase in the critical driving force needed for initiating supercritical crack 
propagation or, equivalently, of the effective fracture energy. The magnitude of the 
increase can be estimated by averaging the fluctuations of the specific fracture energy 
over the length A c as evaluated for the critical reference crack (I = l c ): 

^^{^r^^iff 2 ^- <22) 

At a given applied stress, this increase in the critical driving force implies that the 
critical crack length /* in the presence of disorder will somewhat exceed the critical 
crack length l c in the homogeneous reference system. However, as will be shown below, 
this increase is in general small. 

Above the scale A c , the crack front develops irregular fluctuations. In this regime, 
the depinning theory for manifolds with long-range elasticity predicts the crack shape 
to be self-affine with roughness exponent ( ~ 0.47 [12]. According to this theory, self- 
affine scaling is expected to extend up to a correlation length which diverges as the 
crack driving force approaches the critical value where the crack depins: At this point, 
the crack is expected to be self affine on all scales above A c . However, this is not true 
here, since the nature of the elastic energy functional imposes an intrinsic limit on 
the extension of the self-affine scaling regime. As soon as we consider scales A that are 
larger than the crack length I, Eqn. (JTTj) needs to be replaced by Eq. ( fl8i) for calculating 
the elastic energy of a bulge. Repeating the steps in Eqs. (|T9|) to (I2TI) . we find that 
roughening is restricted to scales below 

A c * « f - (23) 

Above this scale, the elastic interactions again suppress the roughening. Depending on 
the magnitude of the disorder, this creates two possible scenarios (we assume that I x /I y 
and l/l c are of the order of one): 

• For £// > Cy (weak disorder): A c > I > A*; no roughening occurs. 
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• For £/L < Cy (strong disorder): A c < / < A*; the crack roughens on scales between 
A c and A*. 

With roughness exponent ( and upper correlation length A*, we expect the characteristic 
width of a depinning crack to be given by 




Tin) 



(Wr«f ? «e 7 C#. (24) 



4. Comparison with simulation results 

To test the scaling arguments developed in the previous section, we simulate the 
evolution of cracks under conditions of load control. We use a lattice automaton 
technique where we evaluate the displacements at discrete sites (xi,yj) on a two- 
dimensional lattice with lattice constant £. Instead of brittle behavior (t(u) = W{5(u)) 
we consider a triangle-shaped stress-displacement characteristics 

u 

T M , U < U , 

U ° U 

r M (2 ), u <u<2u Q , (25) 

u 

, u > 2uq . 

The fracture energy is here related to the peak strength tm and characteristic 
displacement-to-fracture uq by Wf = tmUq. Disorder is introduced by considering the 
peak strengths (or equivalently the fracture energies) at different sites as independent 
Weibull distributed random variables. In the following we use nondimensional variables 
defined through 

EXT 

r = 7 l T ,E = ^ T ,t/ = ^-,x = f,y=|, (26) 

{tm} {tm} u 4 ? 

such that the scaled average peak strength (T) and fracture energy (W) are by definition 
equal to 1. The scaled crack width is accordingly defined as L — Furthermore, we 
make the simplifying assumption that the mode-II and mode-Ill interaction constants 
are equal, I x — I y — I. The equilibrium equation then reads 



(d 2 U d 2 U\ 

where the non-dimensional interaction constant is J = (/mo)/( t mC 2 )- The non- 
dimensional crack driving force (stress intensity) is given by F = (LS) 2 / (2 J) and 
accordingly the length of the critical reference crack in a homogeneous system is obtained 
from the requirement F c = 1 as L c = \p23 /£. In the simulations we evaluate Uxx an d 
Uyy m terms of the corresponding discrete second order gradients. 

Simulation of a system containing a shear crack is carried out as follows: To create 
a mode-II crack of non-dimensional length 2L, the strength is set to zero over a 'strip' of 
width 2L and U is assumed everywhere zero. Then, the system is loaded by increasing 
the external stress £ from zero in small steps AE. When sites become unstable as 
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Figure 1. System failure stress versus crack length for systems without fluctuations 
and two different values of the interaction parameter J; full lines: analytical 
approximations according to Eqn. (|28[). 



the local (external plus internal) stress exceeds the local strength, the displacement at 
all unstable sites is increased by a small amount AU. Then, new internal stresses are 
re-computed and it is checked again where the sum of the external and internal stresses 
exceeds the local strength. The displacement at the now unstable sites is again increased, 
etc. This is repeated until the system has reached a new stable configuration. Then the 
external stress is increased again and the procedure is repeated until the system fails 
completely {Ui > 2 for all sites). The corresponding critical stress is denoted by Sc. 
The procedure is repeated for different values of AU and AE to ensure that the results 
do not depend on step size. The stress-dependent evolution of the crack front (defined 
as the boundary of the area where C/j > 2) is recorded, as is the failure stress. 

To illustrate how the simulations, which assume semi-brittle behavior (finite values 
of tm and Uq), compare with the analytical results for the brittle limit, we show in Figure 
1 the dependence of failure stress on crack length L for two systems without disorder. 
The failure stresses are well described by the analytical equation 

\/27 

S c = V ._ , (28) 
L + V2J K J 

which for L ^> J approaches the result for the brittle case, E c = a/2J / L. The analytical 
form of Eqn. ( |28i) is in line with the suggestion of Bazant ^20\ to introduce size- dependent 
corrections into Griffith-like criteria for semibrittle materials by replacing the crack 
length with the crack length plus process zone size. In fact, it is easy to see that \fj 
defines the non-dimensional characteristic length of the 'process zone' ahead of the crack 
tip which derives from solving Eq. ( 1271) for the semibrittle law given by Eqn. (l25l) . In 
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comparing simulations and scaling results we shall in the following exclusively consider 
situations close to the brittle limit, L > 20 J, where corrections due to finite process 
zone size are small. 




Figure 2. Crack pinning: Increase in critical stress intensity with increasing disorder; 
initial crack length Lq—40, Lq=80, Lo=160, interaction parameter J = 1; full line: 
AF C oc Cl. 



4-1. Crack pinning 

We now turn to systems with disorder and investigate the disorder-induced increase in 
the critical force required for unstable crack propagation. To this end, we adjust the 
parameters of the Weibull distribution of T such that the average of the distribution 
is kept at (T) = 1 while the coefficient of variation, which serves as a measure of the 
degree of disorder, is varied over a series of simulations. Again we introduce cracks of 
different initial length and increase the external stress on the system until failure occurs. 
This is in general preceded by subcritical crack advances between intermediate pinned 
configurations. We record the mean crack length L = (L(Y)) together with the stress 
such that we can determine the non-dimensional crack driving force (stress intensity) 
as F = ((L(F))S) 2 /(2 J). The disorder-induced increase in critical stress intensity 
at depinning, AF C = F c — 1, is shown in Figure [2] as a function of the coefficient of 
variation of the fracture energy distribution. It is seen that the simulation results are 
in good agreement with Eq. (J22l) which predicts a quadratic relationship between the 
increase in critical stress intensity F c and the coefficient of variation of the fracture 
energy distribution. Also in agreement with Eq. (T221 . it is found that the increase in 
critical stress intensity does not appreciably depend on the length of the crack. 
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4-2. Crack roughening and final crack shape 

We now study how the crack shape evolves during loading. Simulations were carried out 
for cracks of varying initial length Lq in the regime 40 < Lq/ J < 160 where according to 
Figure 1 the brittle approximation can be considered valid. Figure [3] shows a crack shape 
sequence that occurs while the external load is increased until failure. The crack develops 
a complex shape which we analyze in terms of self-affine roughness. To demonstrate 
self-affine roughness we perform a multiscaling analysis where we evaluate the n-th 
order structure function C n (S) := (\L(Y) - L(Y + S)\ n \)^ n (the average runs over all 
positions Y in a given realization and over multiple realizations of the disorder). Self- 
affine behavior is characterized by a power-law relationship C n (S) oc S^ n where ( n = ( 
does not depend on the order of the structure function. 




> 1 1 1 > 1 — 

200 400 600 

Position along crack 



Figure 3. Sequence of crack front positions at increasing external stress; initial 
crack length Lq=80, interaction parameter J = 1, Weibull modulus (3 — 0.5. The 
uppermost curve gives the crack shape immediately before system failure by unstable 
crack propagation. 

Figure H] shows that power-law scaling is indeed observed over about 1.5 orders of 
magnitude in scale. Within the scaling regimes and for 1 < n < 16, we find only a very 
weak n-dependence of the scaling exponent ( n = 0.35 . . . 0.4, justifying the conjecture of 
self-affine behavior. The ( exponent values of about 0.4 are in approximate agreement 
with the theoretical results and experimental observations quoted in the introduction. 
It is important to note that the scaling regimes in our simulations are quite limited: 
Self-affine scaling is limited from below by the pinning length L c and from above by the 
intrinsic correlation length L* = A*/£. 

An analysis of the roughening kinetics has been performed in terms of the lowest- 
order structure factor Cx(S). This is shown in Figure [5] where Ci(S) is plotted for various 
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Figure 4. Multiscaling analysis of crack shapes immediately before failure; parameters 
as in Figure [31 the dotted lines indicate the limits of the scaling regimes used in fitting 
the scaling exponents Q n . 



levels of the crack driving force. As the force increases towards the depinning threshold, 
the crack roughness increases on all scales and a linear scaling regime emerges in the 
double-logarithmic plot. However, it is evident from the figure that the roughening 
kinetics cannot be described by the usual Family- Vicsek form [2T] . 




Figure 5. Evolution of the structure factor with increasing load. Values of the excess 
driving force from bottom to top: AF/AF C = 0.2, 0.4, 0.6, 0.8, 1.0. 
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In spite of this complication, the final shape of the critical crack follows the results 
of the scaling theory outlined in Section 3. Figure [6] shows how the width of critical 
cracks depends on the degree of disorder. In line with Eqn. (Oljl . the width of the 
crack front scales like L 2< *Cy where ( ~ 0.4. The same scaling applies to the difference 
L — Lq of the initial and final crack front positions. This is to be expected since the 
crack roughens by advancing. 




Figure 6. Width of the critical crack (lower symbols, right axis) and crack length 
increase (L) — Lq (upper symbols, left axis) as a function of the coefficient of variation 
of the fracture energy distribution; all lengths have been scaled with L 2 ^ where £ = 0.4; 
the full lines have slope 2( = 0.8; note that the left and right axes have been shifted for 
improved visibility. The error bars represent the variance of values determined from 
10 simulations 



5. Discussion and Conclusions 

We have shown that disorder leads to crack front roughening while simultaneously 
increasing the fracture toughness of the material. Even though the critical stress 
intensity required for unstable crack propagation increases with increasing disorder, 
under conditions of load control this does not necessarily mean that the stress supported 
by the system (with or without crack) increases as well. Figure [7] shows instead a 
complicated behavior: The strength of cracked interfaces first decreases with increasing 
disorder, then increases, then again decreases. To understand this complex behavior 
we have to remember that the presence of disorder allows the crack to advance by 
some amount AL even before it becomes critical. This subcritical advance increases 
the driving force (stress intensity factor) acting on the crack front. At the same time, 
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Figure 7. Critical stress for system failure as a function of the coefficient of variation 
of the fracture energy distribution for different initial crack lengths 



however, the critical force required for crack propagation increases as well by an amount 
AF. Using the condition for the crack to become critical, (SL) 2 /(2J) = (1 + AF), it 
is easy to see that the resulting change AE in critical stress fulfills the approximate 
relation 

AS AF AL 

— < 29 > 

Since AL oc Cy ~ Cy 6 while AF oc Cy, the softening effect of the increasing crack 
length prevails at small disorder, while the strengthening effect of crack front pinning 
takes over at larger disorder, giving rise to the observed non-monotonic behavior. 

However, why does the strength decrease again at very large disorder? This effect 
cannot be understood by considering exclusively the dynamics of the pre-existing crack. 
Rather, one has to allow for the nucleation of new cracks at other sites within the 
system which is facilitated by the disorder as shown in [6J. In the regime of large 
disorder, propagation of the existing crack is inhibited while nucleation of new cracks 
at weak sites elsewhere in the system is facilitated. As a consequence, one ultimately 
reaches a situation where failure occurs not by propagation of the pre-existing crack 
(which instead becomes arrested by the disorder) but by nucleation of one or several 
cracks at other, more favorable sites. 

In conclusion, we have demonstrated that the theory of elastic manifold depinning 
offers a useful framework for understanding the behavior of cracks in disordered 
interfaces underneath thin films. Standard scaling arguments can be used to predict 
the influence of disorder on the critical stress intensity factor required for sustained 
crack propagation, and the geometry of the depinning crack front exhibits features 
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characteristic of a depinning manifold with long-range elastic interactions. However, 
at the same time the standard depinning framework needs to be modified in several 
respects: (i) On large scales, the finite crack length / leads to a change in the elastic 
interactions and to the emergence of an intrinsic correlation length A* > I above which 
the crack remains flat; (ii) The driving force acting on the crack (the stress intensity 
factor) depends on the crack length and can therefore not be controlled independently. 
Subcritical crack roughening proceeds through advance of the crack front which increases 
the stress intensity. As a consequence, one finds a complicated disorder dependence of 
the system failure stress, (iii) At very large disorder, the failure mode changes and crack 
nucleation outwith pre-existing cracks becomes predominant. The system strength in 
this regime is expected to decrease with increasing system size and disorder, similar to 
the behavior of one-dimensional systems studied in [6]. A study of crack nucleation in 
two-dimensional systems remains a task for future investigations. 

Acknowledgments 

Financial support of the European Commission under contract NEST-2005-PATH- 
COM-043386 (NEST pathfinder programme TRIGS) is gratefully acknowledged. 

References 

[1] Alava M J, Nukala P K K N, and Zapperi S, Statistical models of fracture, 2006 Adv. Phys, 55, 
349-476. 

[2] Zaiser M, Slab avalanche release viewed as interface fracture in a random medium, 2004 Ann. 
Glaciol., 38, 79-83. 

[3] Fyffe B, and Zaiser M, The effects of snow variability on slab avalanche release, 2004 Cold Reg. 
Set. Techn., 40, 229-242. 

[4] Heierli J and Zaiser, An analytical model for fracture nucleation in collapsible stratifications, 2006, 

Geophys. Res. Letters, 33, L06501. 
[5] Heierli J, Gumbsch P and Zaiser M, Anticrack nucleation as triggering mechanism for snow slab 

avalanches, 2008, Science, 321, 240-243. 
[6] Zaiser M, Moretti P, Konstantinidis A and Aifantis E C, Nucleation of interfacial shear cracks in 

thin films on disordered substrates, 2009, J. Stat. Mech.: Theory and Experiment, P02047. 
[7] Bouchaud E, The morphology of fracture surfaces, a tool to understand crack propagation in 

complex materials, 2003 Surf. Sci. Review Lett., 10, 797-814. 
[8] Alava M J, Nukala P K V V, and Zapperi S, Morphology of two dimensional fracture surfaces, 

2006, J. Stat. Mech.: Theory and Experiment, 11, L10002. 
[9] Gao H and Rice J R, A first order perturbation analysis on crack trapping by arrays of obstacles, 

1989, J. Appl. Mech., 56, 828-836. 
[10] Ramanathan S, Ertas D, and Fisher D S, Quasistatic crack propagation in heterogeneous media, 

1997, Phys. Rev. Lett., 79, 873-876. 
[11] Ertas D and Kardar M, Critical dynamics of contact line depinning, 1994, Phys. Rev. E, 49, 

R2532-R2535. 

[12] Chauve P, LeDoussal P and Wiese K J, Renormalization of Pinned Elastic Systems: How Does It 

Work Beyond One Loop?, 2001, Phys. Rev. Letters, 86, 1785-89. 
[13] LeDoussal P, Wiese K J, and Chauve P, Two-loop functional renormalization group theory of the 

depinning transition, 2002, Phys. Rev. B, 66, 174201. 



Roughening and pinning of interface cracks in shear delamination of thin films 16 

[14] Delaplace A, Schmittbuhl J, and Maloy K 3, High resolution description of a crack front in a 

heterogeneous Plexiglas block, 1999, Phys. Rev. E, 60, 1337-1343. 
[15] Vcllinga W P, Timmerman R, van Tijum R, and De Hosson J T M, In situ observations of crack 

propagation mechanisms along interfaces between confined polymer layers and glass, 2006, Appl. 

Phys. Lett, 88, 061912. 

[16] Louchet F, Faillettaz J, Daudon D, Bedouin N, Collet E, Lhuissier J and Portal A M, Possible 

deviations from Griffith's criterion in shallow slabs, and consequences on slab avalanche release, 

2002, Natural Hazards and Earth System Sciences, 2 1-5. 
[17] Astrom J A, Alava M J and Timonen J, Roughening of a propagating planar crack front, 2000, 

Phys. Rev. E, 62, 2878-2881. 
[18] Zaiser M, Fyffe B, Moretti P, Konstantinidis A, and Aifantis E C, Pinning and propagation of 

interface cracks in slope failure: ID and 2D considerations, 2004, in: Modelling of Cohesive- 

Frictional Materials, Eds. P.A. Vermeer, W. Ehlers, H.J. Herrmann and E. Ramm, Taylor and 

Francis, London, p. 435-446. 
[19] Arndt P F, and Nattermann T, Criterion for crack formation in disordered materials, 2001 Phys. 

Rev. B, 63, 134204. 

[20] Bazant Z P and Planas J, Fracture and Size Effect in Concrete and Other Quasibrittle Materials, 

1998, CRC Press, Boca Raton (FLA). 
[21] F. Family and T. Vicsck, Scaling of the active zone in the Eden process on percolation networks 

and the ballistic deposition model, 1985, J. Phys. A, 18, L75-L81. 



